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NOMENCLATURE 

(A  partial  list  only) 

h(n)  . 

i.  L. 

the  nl  zero  of  the  Bessel  function  of  order  £. 

f 

the  Coriolis  parameter  2r2sin  cty  . 

k 

the  wave-number  in  the  x-direction. 

P 

the  total  pressure  associated  with  some  particle 
of  water,  except  where  it  occurs  in:  b (n  / , 
or  N  .  In  these  three  cases  it  referspto  p 

the  porder  of  the  Bessel  function. 

p» 

a  small  perturbation  of  pressure. 

t 

time . 

u 

the  velocity  in  the  direction  of  the  x-axis. 

V 

the  velocity  in  the  direction  of  the  y-axis . 

w 

the  velocity  in  the  direction  of  the  z-axis. 

X 

a  horizontal  axis  of  the  Cartesian  coordinate 
system. 

y 

a  horizontal  axis  of  the  Cartesian  coordinate 
system  along  which  there  is  no  change  in  any 
variable . 

z 

the  vertical  axis  of  the  Cartesian  coordinate 
system  aligned  in  the  direction  of  decreasing 
gravity. 

A 

the  ratio  of  the  total  density  variation  change 
to  the  maximum  density. 

JP  ' 

the  Bessel  function  of  the  first  kind,  of  order  £. 

N 

P 

Weber’s  form  of  the  Bessel  function  of  the  second 
kind,  of  order  £. 

c*  __ 

an  inverse  length  scale  associated  with  the  density 
distribution. 

"S  - 

the  vertical  particle  displacement. 

p  - 

the  total  density  of  a  particle. 

V 


NOMENCLATURE 


po 

'Z 

<JO 

n 


the  time-average  density, 
the  time  lag. 

the  latitude,  measured  from  the  equator. 

the  angular  frequency  of  the  free  internal 
gravity  waves. 

the  angular  frequency  of  the  Earth’s  rotation. 
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INTRODUCTION 


1.1  Scope  of  the  work. —  In  recent  years  considerable 
attention  has  been  given  to  the  subject  of  internal  waves. 

Much  of  the  work  has  been  of  a  theoretical  nature,  but  it 
has  become  evident  that  more  data  from  the  open  ocean  is 
necessary  to  provide  a  basis  for  evaluating  the  theory  and 
for  indicating  areas  which  require  further  analytical  treat¬ 
ment.  This  work  is  concerned  with  both  the  theoretical  and 
experimental  aspects  of  internal  waves  in  the  Arctic  Ocean. 

The  theoretical  work  is  an  investigation  of  free  internal 
gravity  waves  in  a  rotating  fluid.  For  a  given  vertical 
profile  of  the  _in  situ  density  the  theory  leads  to  a  descrip¬ 
tion  of  the  corresponding  dispersion  relation  and  eigen¬ 
functions  of  the  vertical  amplitude.  The  mathematical  model 
for  the  density  profile  used  in  this  work  gives  a  good  approx¬ 
imation  to  the  iai  situ  density  profile  found  in  the  Canadian 
Basin  of  the  Arctic  Ocean.  With  suitable  changes  in  the 
parameters  it  could  be  used  to  model  the  density  profile  in 
other  oceans  or  other  fluids. 

The  experimental  work  is  primarily  a  presentation  of  the 
power  spectra  of  internal  wave  records.  These  records  were 
obtained  at  depths  of  60  and  125  meters  during  the  summer  of 
1965.  The  limited  extent  of  the  experiment  precludes  any 
qualitative  comparison  between  the  data  and  the  theory. 
However,  the  experimental  results  should  be  helpful  in  plan¬ 
ning  further  studies  in  the  Arctic  and  in  providing  knowledge 
of  internal  wave  spectra  as  they  occur  in  the  real  ocean. 
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1.2  Brief  1 it erature  survey. —  Stokes(l849)  has  been  credited 
with  the  first  mathematical  investigation  of  waves  on  the 
boundary  between  two  fluids  of  different  density.  A  number 
of  authors  subsequently  elaborated  upon  this  work,  notably 
in  extending  the  treatment  to  a  fluid  with  a  continuously- 
varying  density  profile. 

Work  by  Ekman(l906)  and  Hel land-Hansen  and  Nansen (1909) 
drew  attention  to  the  importance  of  internal  waves  in  the 
ocean.  During  the  early  1900's  investigators  made  further 
attempts  to  learn  about  internal  waves  in  the  ocean,  but  it 
was  Fj eldstad ' s ( 1933 )  analytical  treatment  of  the  subject 
which  provided  the  basis  for  present  work. 

Since  then  there  have  been  numerous  papers,  both  of  a 
theoretical  and  experimental  nature,  attempting  to  explain 
how  internal  waves  are  generated,  how  they  propagate,  and  to 
determine  the  nature  of  the  power  spectrum.  Notable  theoreti¬ 
cal  papers  dealing  with  the  generation  and  propagation  of 
internal  waves  have  been  those  of  Groen(l948),  Haurwitz ( 1948 , 
1950),  and  Rattray ( 1957 , I960) .  Important  papers  of  an  experi¬ 
mental  nature  have  been  published  by  Fjeldstad(l964) ,  Haurwitz 
et  al ( 1959 ) ,  LaFond( 1961 ) ,  Lee(l96l),  and  Reid(l962). 


II.  THEORY  OF  FREE  INTERNAL 
GRAVITY  WAVES  IN  THE  ARCTIC  OCEAN 


2.1  Important  assumptions . —  The  following  assumptions  will 


be  used  in  this  work: 

1.  The  fluid  is  incompressible  and  inviscid. 

2.  The  motions  are  so  small  that  second-order  quantities 
can  be  neglected. 

3.  The  density  variations  are  important  only  where  they 
are  associated  with  the  gravity  term. 


4.  There  is  no  variation  of  any  parameter  in  the  y- 
direction. 


5.  The  density  can  be  separated  into  two  portions,  a 
steady-state  or  average  density  denoted  by  a  zero 
subscript,  and  a  perturbation  density  denoted  by  a 
prime.  That  is: 

-v-  p' 


6.  The  pressure  can  be  separated  into  three  portions, 
a  constant  term  due  to  atmospheric  pressure,  a 
hydrostatic  term  due  to  an  average  density,  and  a 
perturbation  term.  That  is: 

Po  -  ^  po  (V)  i  -v 

7.  The  dynamical  effects  of  the  earth's  sphericity 
can  be  ignored. 

8.  The  only  component  of  the  angular  velocity  of  the 
Earth  which  will  be  considered  is  the  one  opposite 
to  that  of  gravity.  This  assumption  is  valid  only 
for  regions  of  high-latitude  and  mid-latitude. 

2.2  Coordinate  system . —  The  Arctic  Ocean  will  be  approximated 
by  an  ocean  of  constant  depth  D,  on  a  platform  which  is  rotat¬ 
ing  at  a  constant  angular  frequency  t0o.  The  value  of  CG0 

/ 

is  one-half  the  local  Coriolis  parameter  f=2/2sin <P,  where  n 
is  the  Earth's  angular  velocity  of  rotation  and  is  the 
latitude  where  the  experiment  was  performed,  measured  from 
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the  equator.  The  coordinate  system  will  be  the  Cartesian 
coordinate  system  with  the  z-axis  aligned  in  the  direction 
opposite  to  that  of  gravity.  Figure  1  shows  the  coordinate 
system . 

2.3  Equat ions  of  motion . — Following  the  work  of  Fjeldstad 
(1933)  the  above  assumptions  and  coordinate  system  lead  to 
the  following  equations  of  motion: 


<Du 

at 

-  -far 

=  -£•  t?' 

Po  JX 

(2.3.1) 

c)aS 

at 

r  O 

(2.3.2) 

c)axt 

at 

"f. 

(2.3.3) 

au 

•v  a  nS 

4  c)> ur 

■=  o 

(2.3.4) 

¥ 

4- 

yUr  ipo 

r.  O 

(2.3.5) 

In  the  following  work  the  vertical  velocity  w  will  be  re¬ 
placed  by  the  vertical  displacement  %  ,  by  making  use  of  the 
formula : 

dl  =  =>GU-  (2.3.6) 

dt  <Dt 

2.4  The  wave  solutions . —  Since  the  depth  over  the  entire 
basin  has  been  assumed  to  be  uniform  and  the  variation  in  the 
y-direction  has  been  assumed  to  be  zero,  a  solution  to  the 
system  of  equations  (2 . 3 . 1 )- ( 2 . 3 . 5  )  can  be  found  which  has 


the  form: 
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Figure  1„  The  coordinate  system. 
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u  =  U( z ) 
v  =  V(z) 

=  ^)(z) 
p'  =pJT(z) 

p  =  f^2) 


i  (kx  +  oo t ) 

e  v  ' 


(2.4.1) 


Substituting  the  values  of  the  variables  given  in  equation 

(2.4.1)  into  the  equations  ( 2 . 3 . 1 )- (2 . 3 . 5  )  and  eliminating  all 
but  the  vertical  displacement  S  leads  tc  the  following  ordi¬ 
nary  differential  equation: 


X  cL 
po  dt 


=  O 


(2.4.2) 


A  more  convenient  form  of  the  equation  is  obtained  by 
the  substitution: 


which  results  in  the  equation: 


(2.4.3) 


(2.4.4) 


In  the  Arctic  Ocean,  as  in  most  oceans,  the  first  two 
terms  in  the  brackets  of  equation  (2.4.4)  can  be  neglected 
(see  Appendix  I).  A  general  solution  to  (2.4.4)  can  then 
be  obtained  if  the  form  of  the  average  density  distribution 
is  specified.  Numerous  models  for  this  density  distribution 
have  been  investigated.  A  particular  model  which  has  not 
been  investigated  and  gives  a  close  approximation  to  the 
actual  iai  situ  density  distribution  in  the  Arctic  Ocean 
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is  the  following: 


Substituting  this  into  equation  (2.4.4)  gives: 

o 

where  and  ip1-  • 

Equation  (2.4.6)  reduces  to  Bessel’s  equation: 

.1.  ^Znr\  ^  ^ 


d  r 


4-  y  +  (  V^-  ^  ^  W\  =  O 


<v 


(2.4.5) 


(2.4.6) 


(2.4.7) 


by  the  substitution  v^^VYYv.(f|,  where  o<  r  =2^  £  and  —£^_. 

The  general  solution  of  equation  (2.4.7)  is  then: 


yv\  (  y')  —  C-  C  -v  C  Y  ^ 


(2.4.8a) 


or, 

CTt(^e^Vo^(^et)  (2.4.8b) 

2.5  Boundary  conditions . —  To  specify  the  problem  completely, 
boundary  conditions  are  needed.  Most  theoretical  analyses 
of  wave  motion  in  the  deep  ocean  require  that  the  surface 
(z  =  0)  pressure  be  continuous  with  the  atmospheric  pressure. 
However,  the  pack  ice  covering  the  Arctic  Ocean  acts  like  a 
lid,  particularly  for  the  shorter  period  waves,  as  Hunkins 
(1962)  has  shown.  This  means  that  the  dynamic  boundary  condi¬ 
tion  at  the  surface  can  be  replaced  by  a  kinematic  boundary 
condition  requiring  no  vertical  motion  at  the  surface.  For 
the  very  long  waves  in  the  deep  ocean  the  pack  ice  is  not  such 
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an  effective  lid,  but  the  amplitudes  of  the  vertical  motion 
at  the  surface  are  so  small  compared  to  the  wave-length  that 
the  dynamic  boundary  condition  and  the  kinematic  boundary 
condition  are  very  nearly  equivalent. 

The  boundary  condition  at  the  bottom(z  =  -D)  requires 
that  the  particle  displacement  there  be  zero. 

Symbollically,  the  above  arguments  can  be  stated  as: 

3  =  0,  at  z  =  0  (2,5.1) 

*5  =  0,  at  z  =  -D  (2.5.2) 


2.6  The  e igen-values . —  With  the  above  boundary  conditions 
the  eigen-values  are  determined  by  finding  the  roots  of  the 
equation : 

3^ ( (2.6.1) 

'  cx  '  C< 

Writing  equation  (2.6.1)  in  the  form: 


Xf  (V^  r  O  (2.6.2) 


where  ^  -  0  and  k>-  it  can  be  shown  tha  the 


-cxO 

and  h, 

*■»  cx 

roots  of  this  equation  approach  the  roots  of: 


)  -  O  (2.6.3) 

as  >0.  A  demonstration  of  this  is  found  in  Kline(l948), 
who  shows  in  the  same  paper  that  the  roots  of  equation  (2.6.2) 
are  unique  and  positive  as  ^->0,  and  that  the  correspondence 
between  the  roots  of  equations  (2.6.2)  and  (2.6.3)  is  one-to- 


one  . 
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To  determine  the  error  in  replacing  equation  (2.6.2)  by 
(2.6.3)  consider  a  depth  D  of  3000  meters,  which  is  typical 
for  the  Canadian  Basin,  and  a  value  of  CX.  equal  to  0.8xl0~ 


—  1  -5 

meters  ,  then  ^  is  about  10  .  The  most  slowly  converging 


root  of  (2.6.2)  is  the  first  root  of  zero  order  (according  to 
figure  1  (page  34)  of  K1  ine  ( 1948  ).)  According  to  the  same 
reference  the  error  which  results  from  replacing  equation 


(2.6.2)  by  (2.6.3)  for ^  =  10  5  is  of  the  order  of  one  part 
4 

in  10  ,  for  the  first  root  of  the  zero  order.  The  error  for 


all  other  roots  would  be  less  than  this. 

Since  the  maximum  error  for  the  conditions  in  the  Arctic 
Ocean  is  so  small  it  is  reasonable  to  use  equation  (2.6.3)  to 
determine  the  eigen-values.  By  determining  the  values  of  the 
constants  in  equation  (2.4.5)  it  is  possible  to  describe  the 
dispersion  relation  by  making  use  of  equation  (2.6.3).  As  can 
be  seen  from  figure  2,  the  following  values  substituted  into 
equation  (2.4.5)  give  a  close  approximation  to  the  actual  in 
situ  density  distribution  of  the  Arctic  Ocean: 


oc  =  0.8x10 


-2 


meters 


-1 


A  =  5.0x10 


-3 


(2.6.4) 


grams/ cc 


If  b^P)  is  the  n*'*1 
n 


root  of  equation  (2.6.3)  then: 


(2.6.5) 


Substituting  the  necessary  values  of  the  constants  from 
equation  (2.6.4)  into  the  definition  of  given  previously 
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-300- 
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and  making  use  of  the  same  definition  in  equation  (2.6.5) 
leads  to  the  following  dispersion  relation: 

where  60r>an(i  f  are  expressed  in  radians/sec.  and  in 


co 


^  =  -f 


(2.6.6) 


-1 


a 


meters  ,  and,  as  before,  ^ 

Figures  3,4,5,  and  6  show,  respectively,  the  first  three 
modes  of  Sv/Cn  for  the  angular  frequencies  i.4xl0~^  radians/ 

_3  _2 

sec.,  2.0x10  radians/sec.,  and  1.1x10  radians/sec.,  and 


the  dispersion  relation,  equation  (2.6.6). 
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5*/ 
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III.  EXPERIMENTAL  WORK 


3.1  Development  of  the  experiment. —  With  the  exception  of  the 
work  by  Haurwitz  et  al(l959)  measurements  of  internal  waves  in 
the  real  ocean  have  been  obtained  from  stations  located  along 
the  continental  shelf  rather  than  in  the  deep  ocean.  In  an 
attempt  to  add  to  the  information  from  deep  ocean  stations 
an  experiment  was  conducted  on  ice  island  T-3  during  the 
summer  of  1965.  The  ice  island  was  in  the  Canadian  Basin  of 
the  Arctic  Ocean  and  its  position  at  important  times  during 
the  summer  are  given  in  table  1. 

The  Arctic  Ocean  was  chosen  for  the  experiment  for  a 
number  of  reasons.  (l)  The  ice  island  provides  a  very  stable 
platform  from  which  to  conduct  an  investigation  of  waves. 

(2)  The  horizontal  gradients  of  temperature  and  salinity  are 
small.  (3)  There  are  no  swift  currents  which  might  affect 
the  propagation  of  internal  waves.  (5)  The  density  distri¬ 
bution  is  such  that  it  can  be  modeled  quite  well  by  a  simple 
mathematical  expression. 

The  primary  purpose  of  the  experiment  was  to  determine 
the  power  spectrum  of  internal  waves  at  various  depths. 

To  determine  the  amplitude  of  the  internal  waves  the  pro¬ 
cedure  has  been  in  most  cases  to  determine  the  variation  of 
some  depth-dependent  conservative  property  as  a  function  of 
time.  That  is,  if  @(x,y,z,t)  is  a  conservative  property, 
then : 

d®  -  d©  -4.  UL  c)©  r\S  0©  4- mjt  <D©  -  q 
c)*t  <Dk  c)  i 


(3.1.1) 
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Table  1.  Positions  of  T-3. 


DATE 

TIME* 

LATITUDE 

LONGITUDE 

1 7/2 0/6 5 

1200Z 

76°  07.6’  N. 

141°  59’  W. 

2 8/1 0/6 5 

2000Z 

75°  30. 5f  N. 

140°  16 ’  W. 

3 8/2 4/6 5 

2015Z 

75°  57.0’  N. 

139°  38*  W. 

9/  6/65 

2200Z 

75°  24.3’  N. 

139°  51’  W. 

9/10/65 

0130Z 

75°  22. 8»  N. 

140°  10’  W. 

9/10/65 

0800Z 

75°  21.5*  N. 

140°  12 »  W. 

9/11/65 

0600Z 

75°  21.0’  N. 

140°  15’  W. 

9/12/65 

0630Z 

75°  22.2'  N. 

140°  15’  W. 

9/13/65 

0730Z 

75°  23.4’  N. 

140°  14’  W. 

*Date  on  which  hydrographic 

o 

Date  on  which  hydrographic 

3 

Date  on  which  hydrographic 


station 

1-65 

was 

made 

station 

2-65 

was 

made 

station 

3-65 

was 

made 

*Z  indicates  Greenwich  Mean  Time 
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This  property  can  be  considered  as  being  composed  of  a  time- 
average  component  ©)o(x,y,z)  and  a  small  perturbation  compo- 

t 

nent  ©  (x,y,z,t).  If  the  time-average  component  has  large 
characteristic  lengths  in  the  horizontal  then  equation  (3.1.1) 
can  be  written  approximately  as: 


<D©  c)  ©o  =  O  (3.1.2) 

0  "t  c)"t-  ^ 


Integrating  this  equation  leads  to  the  result: 


=  -  d<3>. - 


(3.1.3) 


dh 


from  which  the  amplitudes  for  the  spectral  analysis  can  be 
determined . 

In  the  Arctic  Ocean  one  of  the  variables  which  satisfies 
the  above  requirements  is  the  _in  situ  temperature.  It  is  also 
true  that  the  instrumentation  required  to  measure  in  situ 
temperature  is  much  less  complex,  more  stable  and  more  accu¬ 
rate  than  that  required  to  measure  other  conservative  proper¬ 
ties  such  as  salinity  or  jji  situ  density. 

3.2  Instrumentation. —  The  sensors  used  to  determine  the  ver¬ 
tical  displacements  were  Fenwal  thermistors  (type  GB32P28) 
mounted  in  plastic  tubes  with  a  potting  compound  (Scotch  Cast) 
to  protect  the  joint  where  the  thermistor  was  soldered  to 
the  marine  connector. 

The  thermistor  was  attached  to  the  end  of  a  two-conduct¬ 
or  armored  cable.  A  bight  was  taken  in  the  cable  and  a  40- 
pound  lead  weight  was  attached.  Two  cables  with  thermistors 
attached  in  this  manner  were  lowered  through  a  hole  in  the 
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pack  Ice  near  the  edge  of  the  island.  The  lines  were  measured 
at  the  surface  and  marked  so  that  one  thermistor  would  be  lo¬ 
cated  at  a  depth  of  125  meters  and  the  other  at  a  depth  of 
60  meters.  A  pressure  transducer  was  connected  to  the  end  of 
the  125-meter  cable  to  monitor  depth  variations  caused  by 
changes  in  the  wire  angle.  The  60-meter  cable  was  lashed  to 
the  125-meter  cable  so  that  any  variations  in  the  depth  of 
the  60-meter  cable  could  be  detected  by  the  pressure  trans¬ 
ducer.  The  final  positioning  of  the  equipment  is  shown 
schematically  in  figure  7. 

From  the  hole  in  the  pack  ice  the  leads  went  to  an  equip¬ 
ment  shack  on  the  island,  about  75  meters  from  the  hole. 

Inside  the  shack  was  the  recording  equipment  which  consisted 
of  an  unbalanced  Wheat stone  bridge,  a  Varian  6-10  amplifier 
and  a  large  cylindrical  drum  upon  which  the  recording  paper 
was  placed.  An  ink  pen  driven  by  the  amplifier  recorded  the 
output  on  the  paper  strapped  to  the  drum.  The  drum  revolved 
at  the  rate  of  one  revolution  per  hour  and  translated  at 
about  four  centimeters  per  hour.  At  these  speeds  each  sheet 
of  paper  produced  a  record  approximately  12  hours  in  length. 
Figure  8  is  a  schematic  diagram  of  the  recording  equipment. 

The  motor  driving  the  drum  was  of  the  synchronous  type, 
and  since  the  power  provided  by  the  camp  generator  varied  be¬ 
tween  57  and  62  cycles  per  second  it  was  necessary  to  use  a 
time  mark  provided  by  a  precision  chronometer.  The  chrono¬ 
meter  was  insensitive  to  both  voltage  and  frequency  changes. 

At  the  end  of  each  12-hour  record  the  thermistor  was 


20 


Figure  7.  Schematic  diagram  of  thermistor  array 
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Recording  drum 


Figure  8 


Schematic  diagram  of  recording  equipment 
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removed  from  the  bridge  and  replaced  by  a  high-precision 
decade  resistance  box.  By  varying  the  resistance  of  the 
decade  box  a  calibration  curve  was  obtained  for  each  record. 

The  output  from  the  pressure  transducer  was  recorded  on 
a  Rustrak  recorder  provided  by  the  manufacturer  of  the  trans¬ 
ducer.  The  chart  was  calibrated  in  meters.  Two  scales  were 
provided,  an  absolute  scale  which  gave  the  actual  depth  and 
could  be  read  to  an  accuracy  of  4-2  meters,  and  an  expanded 
scale  which  measured  the  variations  from  some  nominal  depth. 

The  expanded  scale  could  be  read  to  an  accuracy  of  jf.4  meters. 
3.4  Data. —  The  portion  of  the  experiment  with  which  this 
work  is  concerned  covered  the  period  from  0725Z  September  8, 
1965  to  2305Z  September  12,1965.  This  period  was  chosen 
because  the  low  relative  velocity  between  island  and  ocean 
at  the  time  produced  no  observable  wire  angle,  according  to 
pressure  transducer  data.  This  was  also  the  longest  period 
for  which  all  the  equipment  worked  properly. 

It  was  not  possible  to  make  determinations  of  hydrographic 
conditions  during  this  period  due  to  a  malfunction  of  the  sali- 
nometer.  It  was  therefore  necessary  to  assume  that  data 
obtained  from  hydrographic  stations  made  a  few  weeks  earlier 
would  provide  sufficient  knowledge  of  conditions  at  the  time 
of  the  experiment.  This  assumption  was  justified  by  the  fact 
that  the  island  moved  very  little  between  the  time  of  the 
first  hydrographic  station  and  the  end  of  the  experimental 
section  under  consideration  (see  table  l).  It  is  also  well 
known  from  the  work  of  such  investigators  as  Coachman ( 1963 ) 
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and  Worthington  (1953)  that  the  horizontal  gradients  of 
temperature  and  salinity  In  this  region  are  very  small. 

This  is  further  emphasized  by  the  fact  that  the  three  stations 
used  in  this  work  (see  figure  9)  show  very  little  difference 
in  values  of  temperature  and  salinity  as  a  function  of  depth, 
even  though  they  span  a  period  of  several  weeks. 

3.5  Reduction  of  the  data. —  The  resistance  of  the  thermistor 
as  a  function  of  time  was  determined  from  the  record  by  making 
use  of  the  calibration  points  at  the  end  of  each  record 
and  time  marks  provided  by  the  precision  chronometer.  The 
zero  line,  or  null  point,  for  each  record  was  also  obtained 
from  the  time  mark.  This  was  possible  because  the  time 
mark  was  produced  by  shorting  the  input  to  the  G-10  ampli¬ 
fier  (see  figure  8)  at  exactly  one-minute  intervals. 

The  temperature  gradient  with  respect  to  thermistor 

resistance ,  was  obtained  from  the  calibration  curve 
d  R0 

for  each  thermistor,  and  the  temperature  gradient  with 

respect  to  depth,  ,  from  the  hydrographic  stations.  The 

o  1 

particle  displacement  as  a  function  of  time  and  space  was 
then  determined  from  the  equation: 


(3.5.1) 

where  Rf(x,y,z,t)  are  the  perturbations  in  thermistor  resis- 

j  d  > 

tance  obtained  from  the  record,  and  *  =  dT. * — °. 

5ft. 

The  raw  data  was  smoothed  by  eye  such  that  variations 
of  periods  less  than  two  minutes  were  eliminated.  The  records 
could  be  read  to  an  accuracy  of  -1  ohm.  In  terms  of  particle 
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displacement  this  would  correspond  to  +13  centimeters  for  the 
60-meter  record  and  jk52  centimeters  for  the  125-meter  record. 

Figures  10  through  27  are  copies  of  the  original  data. 

3.6  Method  of  analysis . —  Analysis  of  the  data  was  accomplish¬ 
ed  using  methods  described  in  Blackman  and  Tukey(l958). 

For  a  function  of  time  X(t)  generated  by  a  random  process 
the  autocorrelation  function  R  CC)  is  given  by: 


(3.6.1) 


The  Fourier  cosine  transform  of  (3.6.1)  determines  the 
power  spectrum  Px(ca>): 

SCO 

cos(u>/t')  A'X:  (3.6.2) 

o 

When  the  data  points  are  obtained  at  finite  intervals 
the  power  spectrum  is  determined  by  numerical  integration  of 
equations  (3.6.1)  and  (3.6.2).  The  estimate  of  the  autocor¬ 
relation  function  is: 

^  V*vJ!  >  -#.0,1,1,— -S  (3.6.3) 

where  nQ  is  the  total  number  of  data  points  and  s&t  determines 
the  maximum  lag  ^  .  The  raw  estimate  of  the  power  spectrum 
is  given  by: 

P(i£fcV  ^  pt  l.OA-S  (3.6.4) 


P. 

where  £ 


■iv 


o  <  JL  <•  s 

*  “  Vi.  Si  =. 

Since  the  record  is  of  finite  length,  only  estimates  of  smooth¬ 
ed  values  of  the  true  spectral  density  can  be  found. 
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This  is  accomplished  by  smoothing  the  raw  estimates  by  the 
"hamming”  method  described  in  Blackman  and  Tukey(i958): 


lo'  =  P^  +  .40,  P»  0 
5Pxai  =  P»U'V*  .-5H  P.U'  *  P^*  ^ 

s p.l°  *  ‘AC»P»lS'V">*  S4  piV> 


(3.6.5a) 

(3.6.5b) 

(3.6.5c) 


Estimates  of  the  power  spectra  were  calculated  on  an 
IBM  7094  using  a  program  incorporating  the  above  principles. 
The  program  was  made  available  by  M.  Eona  of  the  University 
of  Washington's  Computer  Center. 

Figures  28  through  45  are  estimates  of  the  power  spectrum 
for  each  record. 
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Figure  11.  Record  48c.  2040Z  September  8  -  0800Z  September 

Thermistor  depth  -  60  meters. 
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Figure  13.  Record  50c.  2100Z  September  9  -  0730Z  September  10. 

Thermistor  deDth  -  60  meters. 
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Figure  15.  Record  52c.  2145Z  September  10  -  0850Z  September  11 

Thermistor  depth  -  60  meters. 
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Figure  17.  Record  54c.  2200Z  September  11  -  0900Z  September  12. 

Thermistor  depth  -  60  meters. 
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Figure  20*.  Record  48a.  2040Z  September  8  -  0800Z  September 

Thermistor  depth  -  125  meters. 
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Figure  22.  Record  50a.  2100Z  September  9  -  0730Z  September  10. 

Thermistor  depth  -  125  meters. 
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Figure  24.  Record  52a.  2145Z  September  10  -  0850Z  September  11. 

Thermistor  depth  -  125  meters. 
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Figure  26.  Record  54a.  2200Z  September  11  -  0900Z  SeDtpmbpr  i? 

Thermistor  depth  -  125  meters. 
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Figure  ^>\.  Power  spectrum  record  50c. 
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Figure  ‘^>'2.  Power  spec+rum  op  record  5\ 
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Figure  ‘5*5.  Power  spectrum  record  52_e. 
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Figure  35.  Power  spectrum  o.f  record  5 He. 
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Figure  36.  Power*  spectrum  of  record  55c. 
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Figure  ^40.  Power  spec 


50 


O 

6 


00 

8 


i9 

8 


o 

0 


rJ 

o 

o 


< V 


3 


£ 


d 

in 


51 


d 

(0 

in 

-o 

o 

0 

a) 

v 


0 

£ 

o 

j. 

H- 

0 

<u 

o- 

<P 


0 

0- 


0) 

D 

<r> 

iC 


52 


0 

5 


4) 

0 

0 


id 

o 

o 

« 


0 

o 


cJ 

o 

o 


d 
r n 
in 

— o 

V 

o 

o 

0) 


M- 

0 

8 

D 

v 

U 

<v 

O- 

(p 


0 

Cl 

(0 

T 

4J 

V 

D 

cn 

in 


53 


0 

q 


(0 

0 

o 


* 

8 


(vl 

O 

O 


cf 

in 


V 


0 


o 

Q) 


0 


6 

D 

H- 

o 

<u 

O- 

(0 


<v 

* 

0 

:r 


a) 

V 

0 

c n 

Lu 


54 


-5 


£ 


6 

m 

in 

-n 

o 

o 

(V 

L. 


£ 

0 

V 


15 

0) 

a. 

a) 


* 

o 

fl- 

iP 

T 


i. 

o 

cT> 


Li 


IV.  DISCUSSION  OF  RESULTS 


4.1  Decrease  of  power  with  increasing  frequency. —  An  inter¬ 
esting  feature  of  the  spectral  estimates  (figures  28-45)  is 
the  rapid  decrease  of  power  with  increasing  frequency.  For 
frequencies  greater  than  0.2  c.p.m.  at  60  meters  (figures 
28-36)  this  decrease  is  approximately  proportional  to  f 
For  frequencies  greater  than  0.035  c.p.m.  at  125  meters 

(figures  37-45)  the  decrease  is  somewhat  less,  and  is  approx- 

—5 /2 

imately  proportional  to  f  '  .  This  same  phenomenon  has  been 

observed  by  others,  notably  Haurwitz  et  al(l959).  As  these 
authors  have  pointed  out,  the  actual  decrease  may  be  even 
sharper  for  a  number  of  reasons.  (l)  The  smoothing  process 
(equations  (3 . 6 . 5a , b , c  ) )  used  to  determine  the  power  spectral 
estimates  results  in  a  diffusion  of  power  from  the  high  points 
to  the  low  points.  (2)  Errors  in  picking  off  the  values  from 
the  raw  data  result  in  noise  which  may  increase  the  high  fre¬ 
quency  values  more,  in  proportion,  than  the  low  frequency 
values . 

There  are  a  number  of  possible  explanations  for  this 
decrease  of  power  with  increasing  frequency.  Groen(l948) 
has  shown  that  in  a  rotating,  stratified  fluid  there  is  an 
upper  and  lower  limit  on  frequencies  of  internal  gravity  waves 
which  can  propagate.  The  lower  limit  is  determined  by  the 
value  of  the  local  Coriolis  parameter  f  =  2.QsincP,  the  upper 
limit  by  the  maximum  value  of  the  Brunt-Vdisdia  frequency 
N  The  dispersion  relation  (equation  (2.6.6)) 

Po  dr 
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obtained  previously  is  consistent  with  this,  and,  as  can  he 
seen  from  figure  6,  the  upper  and  lower  limits  on  frequency 
in  the  Arctic  Ocean  are  2.0xl0-2  and  1.4xl0-4  radians  per  sec¬ 
ond,  respectively.  Therefore,  one  would  expect  very  little 

_o 

power  in  frequeniies  greater  than  2.0x10  radians  per  second 
( -19  c.p.m. ). 

The  work  of  Hunkins(l962 )  also  provides  a  basis  for  ex¬ 
plaining,  at  least  qualitatively,  the  decrease  of  power  with 
increasing  frequency.  His  work  was  an  analysis  of  surface 
wave  records  in  the  Arctic,  and  showed  that  the  amplitude  of 
the  surface  waves  increased  with  decrease  in  frequency.  If 
the  internal  waves  are  generated  by  the  surface  waves  the  re¬ 
sulting  spectrum  would  be  expected  to  be  similar  to  the  sur¬ 
face  wave  spectrum.  This  argument  assumes,  of  course,  that 
the  coupling  is  the  same  for  all  wave-lengths,  and  also  ignores 
the  effects  of  resonances.  Actually,  it  would  be  expected 
that  the  efficiency  of  the  coupling  would  be  greater  for  the 
longer  wave-lengths  (lower  frequencies). 

4 . 2  The  spectral  peak  at  .08-. 12  c .p.m. —  The  power  spectral 
extimates  of  the  60-meter  records  all  show  a  distinct  peak  in 
the  region  from  .08  to  .12  c.p.m.  In  the  first  three  records 
(figures  28,  29  and  30)  and,  to  a  lesser  extent,  in  the  others 
this  peak  appears  to  shift  with  time  toward  the  higher  fre¬ 
quencies.  This  may  be  due  to  a  Doppler  effect  caused  by  the 
relative  motion  between  island  and  ocean,  which,  for  these 
frequencies,  would  be  of  the  order  of  10%  for  an  island 
motion  of  10  cm. /sec.  During  the  experiment  the  relative 
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motion  was  of  the  order  of  1-5  cm. /sec.,  as  estimated  from 
pressure  transducer  wire-angle.  The  shifting  of  the  spectral 
peaks  with  time  may  also  he  due  to  the  dispersion  of  the  waves 
travelling  from  a  distant  source. 

Hunkins ( 1962 )  states  that  the  surface  waves  in  this  por¬ 
tion  of  the  spectrum  are  generated  by  local  winds  and,  assum¬ 
ing  that  the  internal  waves  are  coupled  to  the  surface  waves, 
the  internal  waves  probably  have  a  similar  origin. 

4.3  The  very  high  frequency  oscillations . —  As  mentioned  in 
paragraph  3<.5  the  raw  data  was  smoothed  by  eye  to  eliminate 
the  very  high  frequency  oscillations.  These  oscillations, 
from  a  visual  examination  only^  are  found  at  frequencies 
greater  than  approximately  2  c.p.m.  As  can  be  observed  in  a 
number  of  the  raw  data  records  (figures  10  and  11,  for  example) 
the  amplitudes  of  these  oscillations  are  of  the  same  order  as 
the  lower  frequency  motions  for  which  power  spectra  were  calc¬ 
ulated.  At  other  times  these  very  high  frequencies  are  barely 
observable.  The  distinct  gap  in  the  spectrum  and  the  random 
manner  in  which  these  very  high  frequency  oscillations  occur 
would  suggest  that  they  are  of  an  essentially  different  nature 
than  the  lower  frequency  oscillations. 

Electrical  noise  in  the  EMF  driving  the  output  pen  may 
have  been  responsible  for  some  of  these  motions.  However, 
when  the  thermistor  was  replaced  in  the  bridge  by  the  decade 
box  the  very  high  frequency  oscillations  were  not  observed. 

It  is  possible  that  the  motion  of  the  thermistor  cables  through 
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the  Earth’s  magnetic  field  could  have  generated  sufficient 
EMF  to  cause  these  very  high  frequency  oscillations.  It  is 
interesting  to  note  that  data  taken  by  A.  Hansen  of  the  Uni¬ 
versity  of  Washington’s  Department  of  Atmospheric  Sciences 
(personal  communication)  on  a  mechanically-driven  "tide  gauge" 
show  similar  motions  in  the  very  high  frequency  range  mention¬ 
ed  above.  This  "tide  gauge"  measured  the  relative  displace¬ 
ment  of  the  surface  level  with  respect  to  the  island.  This 
provided  a  record  of  surface  motion  analogous  to  the  measure¬ 
ments  in  this  work,  which  were  internal  motions  relative  to 
the  ice  island. 

These  very  high  frequencies  may  also  be  a  result  of  tur¬ 
bulence,  but  until  further  experiments  are  conducted  it  will 
be  difficult  to  be  certain. 

4.4  Recommendations . —  The  results  of  this  work  indicate  the 
presence  of  internal  waves  in  the  deep  Arctic  Ocean.  To  de¬ 
termine  the  amplitudes  of  the  individual  modes,  the  kinetic 
and  potential  energy  as  a  function  of  frequency,  effects  of 
friction,  phase  velocity  and  the  source  of  the  internal  waves 
will  require  instrumentation  in  greater  quantity  than  was 
used  in  this  experiment.  A  triangular  grid  of  three  vertical 
arrays  with  five  to  six  thermistors  in  each  vertical  array 
would  perhaps  be  sufficient. 

At  present,  a  major  problem  is  determining  the  position 
and  velocity  of  the  island  with  sufficient  accuracy.  Before 
any  extended  work  is  begun  it  will  be  necessary  to  solve  this 
problem.  A  number  of  methods  are  available,  one  of  which  is 
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being  used  at  present.  The  group  under  Dr.  L.  K.  Coachman 
of  the  University  of  Washington’s  Department  of  Oceanography 
is  testing  a  drift  meter  which  uses  a  12  kilocycle  "pinger" 
resting  on  the  ocean  bottom.  The  difference  in  travel  time 
of  the  "ping"  to  three  hydrophones  located  on  the  pack  ice 
determines  the  island  motion.  This  device  has  been  used  with 
some  success  on  ARLIS  II,  but,  as  yet,  has  not  given  satis¬ 
factory  results  on  T-3 .  Another  method  of  position  control 
under  consideration  is  satellite  navigation,  and  a  receiver 
may  be  available  for  use  in  the  Arctic  within  a  year  or  two. 

Digital  temperature-recording  equipment  for  the  thermist 
or  arrays  will  also  be  necessary.  The  necessity  of  reducing 
the  data  manually  would  be  eliminated  and  longer  records 
could  be  analyzed  rapidly  and  accurately.  This  type  of  equip 
ment  is  available  commercially  at  present,  and  future  work 
in  internal  waves  in  the  Arctic  will  benefit  by  its  use. 
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APPENDIX  I 


Order  of  magnitude  analysis  for  the  vertical  equation. —  Con¬ 
sider  the  terms: 


I-  Oodtf 
»■ 


in .  C^0  f  A 

Then  if  l/oi  is  the  characteristic  depth  of  the  in  situ 
density  profile  and  A  is  the  ratio  of  the  maximum  density 
change  to  the  maximum  density,  the  conditions  under  which 
terms  I.  and  II.  can  be  neglected  with  respect  to  terra  III 
are : 


or : 


For  the  Arctic  Ocean,  typical  values  of 
indicate  that  the  error  in  neglecting 

3 

of  the  order  of  one  part  in  10  . 


the  above  quantities 
terms  I.  and  II.  is 
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